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ABSTRACT 


M— estimator. one kind of robust estimator, has been used 
in SLR data preprocessing. It has been shown that the 
M — estimator has 50% or more breakdown point. 
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INTRODUCTION 

There are three purposes in preprocessing from a. pass of 
raw satellite range measurements: 

a) Correcting system errors for raw SLR data and forming 
observational files; 

b) Fitting a smoothing function to the range residuals 
from the predicted orbi t, re ject ing noises and outliers 
then estimating accurate for this pass; 

c) Forming QL, FL and NP data files. 

The second term is very impotant for data preprocess, because 
the smoothing function will have effect on quality of NP 
data. 

The smoothing function we used is simply a polynomial in 
time. Generaly, the least squares (LS) estimation is:used to 
solve the parameters of polynominal. But, the LS estimation 
is not a robust estimetion. Sometime, there are a large 
number of noises in raw SLR data, especialy those passes are 
in daylight, the solution of the LS estimation will converge 
to false values. 

In this paper, M-estimator, one kind of robust estimator 
has bee'n used in SLR data preprocessing. It ,has been shown 

that the M-estimators has 50% or more breacdown point e * . 
The breakdown point means that, when the probability of 

noises £ .increases to c* .this method will fail. 


M-ESTIMATOR 

The linear equation is writen: 

y.=X^0 +e. (1) 

i i x 


where 


are observations 
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e is the vector of paramaters to be estimated 

X is the vector of coefficients 

e are random errors. 

The M-estimator, cal led Maximum Likelihood Type Estimator, is 
such an estimator which makes the following objective 
function minimum: 

£F{( y r X^9 )/a}=min. (2) 

i= 1 

e are values estimated for e 

a 2 is variance 

F { • } is an even function 

Different objective functions have different M estimator. In 
this paper we used Hampel estimator, here 
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Equation (2) can be rewriten as: 

2 x i* { (yr x I®)/ a} ®o 

i =1 

where 

* {-}=F’ {•}. 

Then 

L-i * 1 J L i = 1 

and 

(r A /a) / (r^cr) . 

From (3) we have 


(4) 

(5) 

(6) 
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and 


i n 

"2 1 ri 2 

<7 = ) r 

N-p Li i 


( 8 ) 


i = 1 


is number of the paramaters estimated. 


When given the starting values Q q and <j q , we can solve Q by 

(5), (6) and (7). The solution is then iterated between (5) 
and (8) until convergence of the object function. 


|U J+1 -U J | / U J <10 


here 


-3 


U= J F{(y.-X i r e)/a} 

i =1 1 1 


is the times of the iteration. 


PROCEDURES 


The predicted and observed ranges' is R and R at each 

C O 

instant of observation T. After atmospheric correction, 
center of mass correction and delay calibration, we have 
the following range residual equation: 

y.=AR. 

1 1 

=a+bp x +e. . (9) x 

Where a, b are range bias and time bias. 
p is the variability of range. 

Reference show a mathod of caculation which have 50% 
breakdown point. 

a) If totel observation data points are N .which are 
divided into n subgroups equaly and every subgroup inclodes 
four data points, as: 

y y y 

n+l J 2n+1 •'an+i 

y y y 

n+2 2n+2 3n+2 



9-60 


n+k 


2n+k 


3n+k 


2n+n 


3n+n 


n=N/4. 

When noise numbers in raw obsrevation data are less then 
N/2, there must be a subgroup in which contains one noise 
point at most. 

b) For the linear model as (9), we can find the linear 
estimated value of b, for any subgroup k: 

k 


V EVi 

1 = 1 

If b is no-bias, we have: 
k 
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and 


where 
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5 = 2 C (C is a constant to be selected) 
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By solving equations (11), we get-r 
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Thus, the residuats of k-subgroup are 
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(1 = 1. 2, 3, 4) 


kl 


«y — b p 
°kl r kl 


c) For each subgroup, the largest and smallest values of 
r are rejeated. And we can get the initial values a , b 

kO kO 

from remained two data points through follows: 


y =a +b p 

kl kO k<V kj 


(j=l, 2) 


d) Then calcuiation the object function of M~ estimate 
used all observations for every subgroups: 

u k = I F ^r\o~\oK ) (k=I - 2 - • • ■ n > 

i * 1 

where F{-} can be taken from (3), and the initial varlue of a 
can be arbitrarily given, for example 0. 5 'meters. 

e) Select the minimum value from U fc (k=l ,-2, . . . n) . 
Suppose k=m, that is 


U =m i n . 

m 


Then a mQ and b mQ . those are taken from m-subgrop. can be used 

as the initial values a .b . It is sure that the a and b 

o O 0 0 

are taken from ’good’ observation points. 


f) Then we can get 


r i =y r a 0 -Vi 


N 

2 _ 1 

a ) r 

0 N - 2 ^ i 
i a 1 


( i = l , 2 N) 


Because a Q , b Q are obtaied by two data points, they have just 

lower accuracy. From (5) to (8) and iterated until 
covergence, the accurate results a. b can be get as above. 


g) A#ter correcting range bias and time bias, we can get 
a polynomial in time as following: 


A p’ =y -a-bp. 

1 1 1 

=a +a t +a t 2 +a t 3 + 

0 1 i 2 i 3 i 

Useing M — estimator, the parameters of polynomial 

a o’ a i' a 2’ a 3 can s °l ve d- 


CONCLUSION AND DISCUSSION 


Comparing with the LS estimator, M-estimator has its 
advantage as follows: 
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a) It can- be preprocess observation data that contain a 
large amount of noises, for example , a pass for LAGEOS 
in daytime are shown in fig 1, (12/20/1991 8:45 UT) . In this 
pass rate of noise is up to 70%. 

I ...» 

b) At same accuracy, the order of polynomial fitting is 
only 4 using M~est imator , whi le the order is up to 6“8 01 
more with LS estimator. Seeing table 1. 

c) Noise mixed at the parts near the both ends of the 
curve can be detected and deleted. 

Besides, comparison with the method of screen“process in 
and LS estimator, one third time is saved with M— estimator. 
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Table 1. 

Comparison for Two Estimators (Lageos) 


Passes 

Order 

M-est imator 
RMS (cm) Points 

Order 

LS-estimator 

RMS (cm) Points 

Y M D H 

92011011 

4 

5. 8 

28 

8 

5. 9 

29 

92011110 

4 

5. 9 

778 

8 

5. 7 

752 

92011120 

4 

5. 1 

452 

8 

5. 4 

419 

92011317 

4 

6. 0 

216 

4 

6. 6 

212 

92011321 

4 

4. 8 

170 

4 

5. 0 

169 

92011416 

4 

4. 6 

169 

8 

5. 1 

169 

92011419 

4 

5. 9 

94 

8 

5. 9 

86 

92011512 

4 

5. 4 

457 

8 

5. 8 

453 

92022216 

4 

6. 2 

187 

8 

6. 1 

183 

92031116 

4 

5. 6 

326 

8 

6. 1 

326 

92041514 

4 

5. 3 

425 

8 

5. 9 

422 

92041616 

4 

4. 9 

419 

8 

7. 0 

417 

92042011 

6 

5. 6 

60 

8 

5. 3 

56 

92042018 

4 

5. 3 

41 

4 

5. 8 

41 

92042613 

4 

5. 1 

583 

4 

6. 4 

585 

92043015 

4 

4. 1 

77 

8 

5. 3 

83 

92050815 

4 

3. 1 

91 

8 

2. 9 

91 

92052114 

4 

2. 7 

212 

16 

2. 7 

210 

92052312 

4 

2. 6 

170 

12 

2. 5 

170 

92060216 

4 

3. 2 

503 

8 

2. 9 

471 






RESIDUALS OF RANGE (M) 


Fig. 1. 

Residual for A Lageos Pass in Daytime 


STATION ID : 7837 SATELLITE: LAGEOS N= 1506 



ELAPSE TIME (SECOND) 
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Lunar Laser Ranging 




